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mately self-consistent resummation of the thermodynamical potentials of 
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1 Introduction 

At sufficiently high temperature and/or density, quantum chromodynamics 
(QCD) should become accessible by perturbative methods due to asymptotic 
freedom. Thermal perturbation theory is however a surprisingly intricate sub- 
ject, and it was only during the last decade that the necessary methodology has 
been developed, and this development is not yet completely finished. An im- 
portant milestone was the hard-thermal-loop (HTL) resummation programme 
introduced in particular by Braaten and Pisarski which generalized the 
static ring resummations that were known to be required in the calculation of 
thermo- static quantities since long |4^, Q to dynamic quantities such as quasi- 
particle properties as well as production and loss rates. 

There are, however, well-known limitations of a fundamental nature which 
present an impenetrable barrier to perturbation theory at a certain order of 
the coupling (depending on the quantity under consideration), caused by the 
inherently nonperturbative chromo-magnctostatic sector of nonabelian gauge 
theories |6^, ^ ^ . This does not mean that thermal perturbation theory is 
completely futile, though, but rather that at particular points certain nonper- 
turbative input is required in addition ||20|| , for example from lattice calcula- 
tions of the inherently nonperturbative 3-d Yang-Mills theory describing the 
self-interactions of chromo-magnetostatic fields. There are in fact even more 
(less well-known) limitations from coUinear singularities that occur in real-time 
quantities involving external light-like momenta such as the production rate of 
real photons from a QCD plasma ||ll|, |, which require improved resummations 
and/or nonperturbative input. 
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But even in those cases where thermal perturbation theory has not yet run 
into one of those barriers, such as the landmark three-loop calculation of the 
free energy of QCD to order Ci'J'^ by Arnold, Zhai, and others 0, |2^, there 
are severe problems caused by extremely poor convergence and strong renor- 
malizaton scheme dependences which seem to render quantitative predictions 
possible only beyond preposterously high temperatures. While some improve- 
ment seemed possible through tricks like Pade approximations Q , the sad 
conclusion appeared to be that thermal perturbation theory was not applicable 
in QCD at temperatures of practical interest p^ . 

Recently, however, it became clear that the problem of poor convergence is 
not specific to QCD, but arises already in such simple theories as massless scalar 
ip^ theory pO| ] and that alternative resummations can be found that greatly im- 
prove the apparent convergence |^ Q. In these lectures, after an introduction 
to hard thermal loops, I will present an approach to the problem of calcu- 
lating thermodynamical quantities in QCD in an approximately self- consistent 
resummation of hard thermal loops and their next-to-leading order corrections 
[ p^ 16, [T^, 18 that appears to work well down to temperatures a few times 



the deconfinement transition temperature and that suggests that at such tem- 
peratures the still strongly interacting QCD may (at least in certain cases) be 
adequately described by weakly interacting HTL quasiparticles after all. 



2 Thermal field theory 

The Feynman rules of thermal field theory |5^ are most easily formulated 
in the imaginary-time (ITF) or Matsubara formalism. The statistical density 
operator e~^^ is then equivalent to a time evolution operator over an imag- 
inary time interval of length /3, the inverse temperature T~^, and the traces 
in {A) = Tv[e~^^ A\ require (anti-)periodic boundary conditions at the ends of 
the imaginary time interval for bosons (fermions). This gives rise to discrete 
imaginary (Matsubara) frequencies when going to momentum space, so that the 
only change in the Feynman rules is in the replacement 

/ ^ - /^-^ E / - m-f5..m. (1) 

with 

, „ .„ f Z bosons 

k^-^^^'Tv^ Z-i fermions ' 

Since it is usually hard to evaluate the resulting sums directly, they are best 
turned into integrals again by writing, in the bosonic case, 

T f{ko = 27rtnT) = ^ / dfco/(fco)f coth^ 
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^ V ' v ' 

vacuum contribution thermal contribution 

for a function /(fco) that is regular for ka G iM and where C is a contour encircling 
only the poles of coth at the location of the Matsubara frequencies. (In the 
fermionic case, an analogous formula can be easily obtained by using tanh in 
place of coth, which in the right-hand side leads to minus the Fermi-Dirac 
distribution in place of the Bose-Einstein one.) 

This evaluation through contour integrals leads to a nice separation into a 
(Wick rotated) T = contribution, and a purely thermal one, which vanishes 
for P —f oo, i.e. T — > 0. There exists also a formulation directly in real time 
(Schwinger-Keldysh and variants thereof) |Q where this separation is already 
conspicuous in the Feynman rules. However, this requires a doubling of fields, 
a 2 X 2 matrix structure for propagators and even more components for n-point 
vertex functions, which in fact correspond to the many possibilities of analytic 
continuation to real frequencies (in particular if there are several independent 
external frequencies) |^ . 

A simple case that leads to a thermal contribution in (^ is the frequently 
occurring one that there is a simple pole in /(/cq) at fco — ±E with E > 0. Then 
one has 

- — |— — = ±{—iyna(E) + vacuum contributions, (4) 
^-^ kn ± E 

n 

where na{E) — [e^^ — a\~^ and a — +1 for bosons, —1 for fermions. 



3 Hard thermal loops 

The simplest example of a hard thermal loop is given by the one-loop self- 
energy diagram in a scalar gLp"^ theory. This is a tadpole diagram, independent 
of external momenta: 

K 



where we have introduced the notation Y,^, =TY.n! -^^^ = ('^"^ '^") 



With E = V k^ + m? , its thermal contribution is easily evaluated as 
-1 _ ^ 1 / 1 1 



^ K^-m^ ^ 2E \ko - E ko + E 

(The vacuum contribution is removed by standard (T^O) renormalization.) 
Without the Bose-Einstein factor, this integral would be quadratically diver- 
gent. Thanks to the former, it is finite, but dominated by momenta k ^ T. If 
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T ^ m, the mass terms of the T=0 theory can be neglected, and the leading 
contribution to the self-energy is given by a self-energy contribution proportional 
to T^, which is called a hard thermal loop (HTL): 

^='-^T.^-s'^ Tdkk n{k) = g-T- , (7) 

K •'^ 

where the hat is a reminder of the HTL approximation - in this case the neglect 
of the bare mass m. 

n clearly corresponds to a mass term for scalar excitations generated by 
interactions with the particles in the heat bath. It should be understood, how- 
ever, that this thermal mass is qualitatively different from ordinary masses. In 
particular, one can show that unlike ordinary masses it does not spoil conformal 
invariance (if any) as it does not contribute to the trace of the energy momentum 
tensor ||63| . 

This example for a HTL is, however, deceptively simple. In general, thermal 
masses are not constant but depend on momentum, that is, they correspond to 
nonlocal terms in an effective action. 

The thermal masses generated for gauge fields are of this form. Indeed, a 
constant mass term would violate gauge invariance. Let us consider as a simple 
gauge theory example the case of scalar electrodynamics with covariant gauge 
fixing, 

C - iD,,4>rD>^4> - \f^uF^^'' - ^{d^A^^f (8) 

where = 9^ -f ieA^. The tree- level propagator reads G"^, = ^^^/iir^ -|- 
(a — 1)K I and this will receive thermal corrections through the photon 
polarization tensor, which in momentum space reads 

H^.(X) = e2^ 
p 

The last term is from the seagull diagram, which is essentially the same as the 
tadpole diagram before. The first term is however dependent on the external 
momentum, and there is no reason to expect a Lorentz invariant form, because 
the heat bath is singling out a preferred frame of reference. 

There are in fact four symmetric tensors that can be built from the available 
quantities 5^1/, isT^, and the four- velocity of the heat bath which we have tacitly 
chosen as = (^JJ- 

In electrodynamics, the polarization tensor has to be transverse, K^H^n = 0. 
This additional requirement still leaves two possible tensors, for from U and K 
one can build the transverse vector [/^ = K^Uf^ — [U ■ K)K^. We choose them 
as 

A^,^g^,-K^K,IK^-B^,, Bf,,^U^U,/U\ (10) 

A^i, can easily be shown to have vanishing components Aq^, = = ^^Oi whereas 
Aij = —Sij + kikj/k'^, so this is a projection onto spatially transverse momenta; 
Bfii, is a projector orthogonal to A^^. 



{2P - K)^{2P ~ K), 2g^, 



P^{P-Kf 



p2 



(9) 
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In nonabelian gauge theories like QCD, one generally has a more complicated 
structure. Then one needs two more, nontransverse, tensors, which one may 
choose as 



Decomposing 11^,^ = -UfA^^ - II^B^,^ - 11^(7^^ - n^Z?^^ and 



one has 



At 
A. 



a 



- an,! 



- aUd 



(11) 

(12) 

(13) 
(14) 



In view of the explicit gauge parameter dependences (there are in fact more 
hidden within the structure functions of 11) , it is remarkable that one can prove 
that the singularities of At and Ai are gauge- fixing independent ^ [z^ . 

In electrodynamics this situation is much simpler (unless one introduces 
nonlinear gauge fixing) : 11^^ is both transverse and completely gauge parameter 
independent. 

Because of transversality, which is easily verified for (^, there are only two 
independent components, e.g. H'^ and IIoo. The former reads 



-4 

p2 - p2(p_^)2 



p. .0(T)^^ 



0{T) 
(15) 

where the term proportional to does not constitute a hard thermal loop, 
because its integrand does not involve a quadratic divergence in its vacuum 
piece. For fco, k^T, W^iK) ~ fl^^ = e^rVS. Notice, however, that this 'HTL 
approximation' remains valid even for ko,k ^ T as long as ^ . 

The other component, IIoo is more complicated to evaluate. Its HTL piece, 
which is contained in 



noo(if) = e- 



E 



4po(po - fco) 



2 

P2 



can however be extracted rather easily using 
1 1.11 



P2 2^po 



P Po+P 



1 



1 



Po + Xpo + Y X-Y'p„ + Y 



Po 



X 



(16) 



which gives a number of terms of the form jyiy ("■(^) ~ '^(^)) after summing 
over the Matsubara frequencies. Now, HTL contributions arise from p T, and 
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for fco, fc <C T one can approximate the energies X and Y by ±p, except when 
two hard energies form a soft difference. This gives 

^ ' n(\p~k\) + n{p)]\ (17) 



kQ-p+\p-k 



where the tadpole-Uke last term in has canceUed against terms where the 
energy denominators contain the sum of two hard energies. 

Since k <^ p ^ T , we have n(\p—k\) — n{j)) ^ n' {p)^p—k\ —p] and |p— fc| —p ~ 
p ■ k/p = zk, so the HTL piece of IIoo is finally given by 



UooiK) ^ 2e2 / ^n'ip) C dz 



fco , fco + fc 
1 m 



fco 



fco — zk 



2k ko — k 



(18) 



Originally, fco was restricted to Matsubara frequencies. In order to allow for 
soft ko ^ T without being restricted to the zero mode, we in fact need analytic 
continuation, e.g. ko lu + ie for retarded boundary conditions, and this defines 
the cut of the logarithm in (p^. 

The above results for 11'^ and Iloo{K) are actually universal. They have 
the same form in nonabelian gauge theories, only the overall coefficient differs. 
In SU(iV) gauge theories with Nf quark flavors one just needs to replace 
[79[ g'^iN + Nf/2). They also retain their form in the presence of a 

nonvanishing chemical potential /i^, which leads to the replacement of —^ 

+ 3/i^/7r^ in the contributions cx Nf. In the HTL approximation, there are 
moreover no nontransverse contributiong|^ so the HTL gauge boson propagator 
involves two independent branches determined by lie = —^looK^/k^ and lit = 

(H';, - n,)/2. 

The poles of the propagators Af and determine the dispersion laws of 
two sorts of quasiparticles, which in contrast to the scalar ip"^ example are not 
given by simple mass hyperboloids. These are displayed in Fig. ^ in a plot of 
ojt.i{k) in quadratic scales (where a relativistic mass hyperboloid would show 
up as a straight line parallel to the light-cone). 

Above a common plasma frequency Wpi. = eT/3, there are propagating 
modes, which for large momenta in the transverse branch tend to a mass hyper- 
boloid with asymptotic mass = i^pi ^ ^"^^ branch i approach the light- 
cone exponentially with exponentially vanishing residue. Indeed, this mode does 
not have an analogue in the T—Q theory but is a purely collective phenomenon, 
so it has to disappear from the spectrum as fc — > oo. The spatially transverse 
mode, on the other hand, represents quasiparticles that are in-medium versions 
of the physical polarisations of gauge bosons. 

^At order kT one has however lie ^ for a ^ 1 in the nonabehan case, and 11^ ^ at 
two-loop order whenever lie 7^ 0. 
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- 



screening | propagation 



Figure 1: The location of the zeros of A^"'"'^^^ (spatiaUy transverse gauge 
bosons) and of A^'^'^^^ (longitudinal plasmons) in quadratic scales such as 
to show propagating modes and screening phenomena on one plot. 



For uj < Wpi., |k| is the inverse screening length, which in the static limit 
vanishes for mode t (absence of magnetostatic screening), but reaches the Debye 
mass, rnjj — Sojpj , for mode I (electrostatic screening). A vanishing magnetic 
screening mass is required by gauge invariance in abelian gauge theories [ p4| , Q , 
but not in the nonabelian case. In fact, lattice simulations of gauge fixed prop- 
agators in nonabelian theories do find a screening behaviour in the transverse 
sector, however the corresponding singularity is certainly quite different from a 
simple pole |ga . 

For uj"^ <k^, there is a large imaginary part ^ e^T^ from ( p^ ) which prevents 
the appearance of poles in this region. This imaginary part corresponds to 
the possibility of Landau damping, which is the transfer of energy from soft 
fields to hard plasma constituents moving in phase with the field j6^, |l3| and 
is an important part of the spectral density of HTL propagators. At higher, 
subleading orders of perturbation theory, it is, however, not protected against 
gauge dependences in nonabelian gauge theories. 

To complete the discussion of HTL's in scalar electrodynamics, let us also 
consider briefiy the scalar self-energy. In contrast to theory, this is now a 
nonlocal quantity. Nevertheless, the HTL part is still a constant thermal mass 
as given by the first term in 



2(a-l)K-P 



(19) 



The other terms are not proportional to because the integrands do not grow 
sufficiently at large momenta. They are even gauge parameter dependent, in 
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contrast to the HTL piece. Notice that these gauge dependent terms vanish on 
the lowest-order mass shell K'^ = 0. 

In QCD (and already in spinor QED), we also need to consider the fermion 
self-energies Q. In the ultrarelativistic high-temperature limit, bare masses 
can be neglected. In this case the fermion self-energy can be parametrized by 

Y.{K) = a{K)-i^ + h{K)k-^/k . (20) 

Using the same methods as above one easily computes the HTL contributions 



1 e'^T'^ 

-tr = uja + kb ^ — ^ = M^, (21) 

1 n e^T^ ^ uj + k 

where in nonabelian gauge theories now g^i^"^ — l)/(2iV). 
The structure of the HTL fermion propagator is 

S{K) = I 7" + 1 A+ + I 7" - ^U- (23) 



with Ajd^ = — [lj =F (fc + S±)] and S± =b±a. The two branches correspond to 
spinors whose chirality is equal (-|-) or opposite (— ) to their helicity. 

The additional collective modes of branch (— ) ( "plasminos" ) have a curious 
minimum of uj at lo/M « 0.93 and \k\/M « 0.41 and approach the light- 
cone for large momenta, but with exponentially vanishing residue. The regular 
branch approaches a mass hyperboloid (in Fig. || a straight line parallel to the 
diagonal) with asymptotic mass \/2M. Again, for space-like momenta, < 0, 
there is a large imaginary part corresponding to Landau damping, which now 
corresponds to the transmutation of soft fermionic fields together with hard 
fermionic (bosonic) plasma constituents into hard bosonic (fermionic) ones. 




4 Standard HTL perturbation theory 

When writing down dressed propagators, as we have done above in the discus- 
sion of the spectrum of thermal quasiparticles, we have already performed a 
resummation of infinitely many loops — the geometric series of self-energy inser- 
tions according to Dyson's equation, which in the scalar theory is simply 

Clearly, the perturbative version is useful only when ^ (gT)^. In particular, 
it fails for fco,fc ^ gT, the scale where collective phenomena transform the 
familiar quanta into quasiparticles. 
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Figure 2: The location of the zeros of A^^ in the HTL approximation in 
quadratic scales. 



The appearance of thermal masses presents unavoidable problems at some 
order of perturbation theory, namely when these higher orders probe the en- 
ergy/momentum scale /co, fc ~ gT, and, in particular in the bosonic sector, the 
Bose-Einstein distribution factors enhance the sensitivity to the infrared. 

It is therefore mandatory to switch from bare perturbation theory to one 
that uses resummed (dressed) propagators. In simple cases like scalar </5^ theory, 
where the only HTL is a constant mass term, this resummation has been studied 
already long ago Conceptually (although not practically) equally simple is 
the example of scalar electrodynamics, where the only HTL's are the self energy 
diagrams considered above. However, already in spinor QED and to a larger 
extent in QCD, it turns out that there are HTL vertex functions which 
have to be treated on a par with the HTL self energies to achieve a systematic 
resummed perturbation theory psf . Indeed, if A^-point vertex functions give 
rise to HTL's cx T^, they are as important as bare vertices when the momentum 
scale is ~ gT: 



N 
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N-2i4-N 



dA 



N 



k^gT 



(25) 



In fact, even for N so high that there is no comparable tree-level vertex, such 
vertex functions have to be resummed. Already in spinor QED, there exist 
one-loop HTL vertex functions with two external fermion lines and an arbitrary 
number of gauge boson lines. In QCD, there are even more HTL vertex func- 
tions without tree-level analogue. It turns out that one-loop diagrams with an 
arbitrary number of external gauge boson lines are HTL. 
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4.1 HTL effective actions 



The HTL resummation programme can be understood as the transition from 
the fundamental Lagrangian to an effective one generated by 'integrating out' 
the hard momentum modes kg, k T in one- loop order. 

In scalar (p'^ theory, this effective theory differs from the bare one only in a 
simple thermal mass term, ^C^^iar = 

For gauge bosons and fermions, the effective theory is necessarily a non- 
local one as gauge invariance forbids a simple local thermal mass term and 
in the case of fermions it is due to the fact that HTL's do not spoil chiral 
symmetries. Remarkably, the infinitely many non-local vertex functions can be 
summarized by a comparatively simple and manifestly gauge-invariant integral 
representation 0, ||, |^ 

J Att ^ V- D{A) ^ 



S f dfly^^^ Vo,V^ 

47r {vDadj.iA)) 



V — {l,v) is a light-like 4-vector, i.e. with v'^ — 1, and its spatial components 
are averaged over by dfl-^. v is the remnant of the hard plasma constituents' 
momenta ~ Tv^, namely their light- like 4- velocity, and the overall scale T 
has combined with the coupling constant to form the scale of thermal masses, 
M, LUpi, gT. 

The covariant derivatives in the denominators of ( p6| ) are responsible for the 
fact that there are infinitely many HTL's involving external fermions and an 
arbitrary number of gauge bosons, even in QED, where only the pure gauge- 
field sector becomes bilinear because of fadjX^) ^ d. 

Technically, resummed perturbation theory amounts to the replacement 

Cel. Cd. + C^^^ - ££^^^ (27) 

where £ is a loop counting parameter that is sent to 1 in the end, after the last 
term has been treated as a 'thermal counterterm'. 

Because C^'^^ has been derived under the assumption of soft external mo- 
menta, this prescription is in fact only to be followed for soft propagators and 
vertices [ p3[ . Those involving hard momenta (if present) do not require this 
resummation, and they can be excluded from this resummation by the intro- 
duction of some intermediate scale A with gT <C A ^ T. Complete results have 
to come out independent of A, of course. 



4.2 Example: NLO terms in scalar electrodynamics 

Massless scalar electrodynamics |Q is a particularly simple toy model as its 
HTL effective action (|2^) is bilinear in all fields. The scalar HTL self-energy 
( p^ is moreover a simple mass term, and in order to consider the one-loop 
corrections to the photon polarization tensor, nothing more is needed. Let us 



10 



consider two limiting cases of this to illustrate some important points of the 
HTL resummation programme. 



4.2.1 Debye mass 



The Debye mass, i.e. the inverse screening length of electrostatic fields, is 
determined by the zero of A^(fco = 0, fc) at imaginary k leading to = 
noo(fco = 0,k^ = — m|,). Because at leading order IIoo(fco = 0,k) turns out to 
be independent of k (see (|l8|)), the frequently found definition ji^ of m|, as 
noo(fco = 0, fc — !■ 0) happens to be correct, but becomes unphysical in general 
[ [71] : beyond LO, only the former, self-consistent definition is renormalization- 
group invariant and (in nonabelian theories) gauge invariant. 

The resummation programme sketched above makes the LO result rh^ = 
part of the new lowest-order Lagrangian £ci. + C^"^^- The NLO cor- 
rection therefore is given by one-loop diagrams using this Lagrangian. These 
are scalar loops which now have massive propagators with thermal mass /i^ = 
S = from (|l9|), from which the HTL result rnjy has to be subtracted as 

thermal counter-term. Indeed, without this subtraction, the LO result would 
be generated a second time, since for large loop momenta the thermal mass of 
the scalar is negligible. Because of the subtraction, the one-loop integrals are 
now receiving their leading contributions from soft loop momenta k ^ eT: 



STloo{0,q) 



dkk^ 



K 



fi 



1 



kq 



In 



2k 



2k 



(m ^ 0) (28) 



To obtain the leading contribution, we can replace n(E) — )■ T/E which gives 



5noo(0,g) = --^y dk I 



fc2 



1 - - In 

q 



2k + q 



2k - q 



+ 0{e^q^ ln(r)). 



(29) 

The NLO correction to the Debye mass is now given by evaluating this at 
q = iiTiD- Incidentially, the above integral is g- independent, as can be seen from 
an integration by parts, which finally gives 



TOn = riin - 



2n 



(30) 



Notice that the perturbative result at NLO involves a single power of e and so 
is non-analytic in a = e2/(47r). 

The above calculation can actually be simplified by noting ||] that the only 
terms capable of producing odd powers in e are the n — terms in the sums 
over Matsubara frequencies in (|2^). For n 7^ 0, one has —k^ = (2TTn)'^T^ /^^ 
so that the thermal masses can be expanded out, leading to powers of only. 
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Keeping only the n = contributions, the first integral in ( psD vanishes, and we 
have 

moo(0,9) = 2e Ttj —————— ^ _ (31) 



(27r)3 



27r 



in agreement with (|30|). Here we have introduced dimensional regularization 
(with mass scale a) to render the integral finite as in Ref. [^. 

As we shall see presently, this simplified resummation by dimensional reduc- 
tion is only possible in the static case (go = 0), and not for dynamical quantities. 



4.2.2 Plasma frequency 

Propagating modes exist for frequencies qa > LUp\, , and this plasma frequency is 
the same for the transverse and longitudinal modes, because it corresponds to 
the long-wavelength limit g — > 0, and with q = there is no way to distinguish 



/3 = e^TV9. 



the polarizations. In the HTL approximation. 

The NLO correction in the case of scalar electrodynamics can be calculated 
in full analogy to (p8|), but now with go = <^p;. and q ^ 0. Because of g = 0, 
the angular integrals are now trivial and one finds 



(5n,(go,0) -(5nt(go,0) 



e^T 



+ ^([^^^-90^/^-^^)}. (32) 



Evaluated at go ~ ^pi. this gives (Jcjpj 



L^pi. mis gives oWpj /cjpj = — e(8\/2 — 9)/(27r) « — 0.37e. 
Notice that without resumming /i, the result would have been completely mis- 
leading: evaluating the unresummed result, i.e. ( ^2|) with /i = 0, at go = Wpi. 
would have given a purely imaginary result that one would have wrongly iden- 
tified with a damping constant .0 

Furthermore, the correct result ( ^2|) is now only obtained if the nonstatic 
modes are resummed along with the static ones. If one keeps only the zero modes 
and ignors that in the imaginary time formalism the external frequency go has 
to be a multiple of ^-kiT (and so cannot be soft and nonzero) , but immediately 
'continues' to go = Wpi. ~ eT, one would find 



^n^,t(9o,o)|o-mod 




90 



(9o - M^) + 



(33) 



which clearly differs from (^^. The resulting Jwpj would in fact be only about 
a quarter of the true result. So dynamic quantities require the full HTL re- 
summation method; resumming only the zero modes is not sufficient (see also 

m 

^In scalar electrodynamics, unlike QCD, the plasmon damping is a higher order effect 
because the scalar HTL quasiparticles are heavier than plasmons and moreover do not have 
Landau damping cuts in the HTL approximation. 
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4.3 NLO corrections for QCD quasiparticles 

The calculation of NLO corrections to the long-wavelength plasmons in QCD 
was in fact one of the first applications of the HTL resummation programme. 
In particular, the damping constant of order g^T ~ gWpi. was the subject of a 
long controversy (in particular with regard to its gauge-fixing (in)dependence) 
before it was calculated in Ref. with the gauge independent]^ result 

7 « 0.264%/iV5Wpi. (34) 

(for pure-glue QCD). The analogous calculation for the damping constant of 
long- wavelength fermionic quasiparticles was carried out in Refs. [|o[ The 
significance of these results is that gluonic and fermionic quasiparticles, which 
in the HTL approximation appear to be stable, experience damping. For g ^ 1, 
they are weakly damped, whereas for g 1 which is more relevant for experi- 
mentally accessible quark-gluon plasmas, the damping is significant: 7 ~ ^^pi.- 
The NLO correction to the gluonic plasma frequency has also been calculated 
with the resuh (Jw^j /w^j « -O.IS/ZV g. 

The NLO correction to the Debye mass, however, runs into IR problems. 
Naively one would expect problems from the masslessness of magnetostatic 
gluons only at two-loop order resummed perturbation theory, when their self- 
interactions become relevant. However, because gauge independence requires 
evaluation on mass-shell (which in the case of the Debye mass means 50 = 0, 
(f = —mjj), there appear 'mass-shell singularities' caused by the massless mag- 
netostatic modes. Because these singularities are only logarithmic, the leading 
log is perturbatively calculable and reads IT^, ItI] 



^ml/rnl^l^^^^glo,l+Oig). (35) 

For small coupling, the logarithm dominates over the non-perturbative constant 
behind the logarithm, and thus the perturbative prediction is that of a positive 
correction to the screening mass. Indeed, lattice simulations of both (gauge- 
fixed) chromo-clcctrostatic propagators J2q] and gauge- invariant lattice defini- 
tions of the nonabelian Debye mass |55|] give significant positive corrections 
to the HTL value.| 

Such a logarithmic sensitivity to the nonperturbative physics of the chromo- 
magnetostatic sector has in fact been encountered earlier in the damping of hard 

•^In fact, later investigations using covariant gauges encountered again gauge dependences 
1^, which are removed, however, when an infrared cut-off is retained while taking the 
on-mass-shell limit |rc| ] . This means that there is an infrared singularity in the residue of the 
pole, which is gauge dependent, but the position of the singulariiy inp longer a simple pole) 
is gauge independent in accordance with the theorems of Refs. [Bll, p2j. 

*These positive corrections are so large in fact that the nonperturbative contributions at or- 
der g dominate over the HTL result in the gauge-invariant lattice definitions. When extracted 
from lattice propagators, which also give gauge-independent results [M, these nonperturba- 
tive contributions are considerably smaller and about 1/3 of them can be accounted for by 
one-loop resummed perturbation theor y if one introduces a simple phenomenological magnetic 
screening mass taken from the lattice |72| . 
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excitations ^ More generally, it arises for all propagating modes |^ 
as well as for all finite screening lengths |Q. For nonzero wave- vector q, only 
the real corrections to the dispersion law u — uj{q) turn out to be IR safe in 
one-loop resummed perturbation theory. 

5 HTL-resummed thermodynamics 

In the previous section we have seen that dynamic quantities cannot be treated 
by the simplified resummation scheme that resums only static modes. For 
static quantities like the thermodynamic potential, however, such a resumma- 
tion works in the sense that it gives a scheme to systematically compute the 
series expansion in powers (and log's) of the coupling. This calculation has 
been performed to order a^^^ in QCD with the result (for pure glue) 

^ = -^/^ = 1 - ^ + 30 - + ^{-f log - 36 

45 [ 47r TT 2 TT TT 

- 1^ (-^ (-4 

Unfortunately, this is very poorly convergent: only when < 0.05 one has 
apparent convergence, but this corresponds to temperatures higher than lO^Tc 
I 

In what follows we shall attempt a different route that resums also the non- 
static modes, and tries to keep resummation effects even when they are formally 
of higher order than that achievable at a given loop order. 

5.1 Screened perturbation theory 

In scalar (p'^ theory, it has been shown Q that the convergence of thermal 
perturbation theory can be improved if the thermal mass of the scalar quasi- 
particles is kept within thermal integrals and not treated as proportional to a 
coupling constant when setting up the perturbation series. Technically, this is 
just as in (p7[), but without the requirement that the resummation has to take 
place in soft quantities only. Because this changes the UV structure at any finite 
order of perturbation theory, this introduces new UV divergences and associated 
renormalization scheme dependences, which in principle can become arbitrarily 
large. But starting from two-loop order, these can be minimized if the thermal 
mass used in screened perturbation theory is determined by a variational prin- 
ciple, i.e. a prinicple of 'minimal sensitivity' to the mass parameter used in this 
reorganization of perturbation theory ^, . 

In Refs. |, |, this approach has been adapted to a one-loop calculation of 
the thermodynamic potential of QCD where in place of a simple mass term the 
gauge-invariant HTL effective action is used. While the leading-order interaction 
term cx is incomplete (in fact, it is over-included), it does contain the plasmon 
term oc without leading to the disastrous result of a thermodynamic pressure 
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in excess of the Stefan-Boltzmann limit. However, it remains to be seen if this 
is stih true in the (technically very difficult) two-loop approximation, which 
presumably has to contribute with positive sign to make up for the over-included 
negative leading-order interaction term of the one-loop approximation. 

In the following, I shall present instead the result of a different strategy. 
Instead of using the HTL effective action as a gauge invariant mass term for an 
optimization of perturbation theory, which does not (have to) care whether the 
HTL effective action remains accurate for hard momenta (which it does not), the 
formalism of self-consistent "^-derivable" approximations will be invoked to 
find expressions that keep resummation effects for both soft and hard momenta. 
As we shall see, the leading-order effects will arise exclusively from kinematical 
regimes where the HTL approximation remains justifiable. Moreover, it will be 
possible to avoid the spurious UV problems of screened perturbation theory. 



5.2 Approximately self-consistent resummations 

In the Luttinger-Ward representation of the thermodynamic potential H, — —PV 
[ |62| (to particle physicists often more familiar as the composite operator effective 
potential [^7| ) is expressed as a functional of full propagators D and two-particle 
irreducible diagrams. Considering for simplicity a scalar field theory for the 
moment, n[D] has the form 



n[D] = -T log Z ^ Trlog D-^ - ^TTtUD + T<P[D] 

= J -0^niu;)lm [log D^^lj, k) ~ U{cj,k)Diu;,k)] +T^D], (37) 

where Tr denotes the trace in configuration space, and $[-D] is the sum of the 
2-particle-irreducible "skeleton" diagrams 

1/12^^ +1/8 (^^Q + l/48^^+... (38) 

The self energy II — D^^ — D^^ , where Do is the bare propagator, is related to 
by 

6<I>[D]/6D = iH. (39) 

An important property of the functional r2[-D], which is easily verified using 
(^9|), is that it is stationary under variations of D: 

6n[D]/6D = 0. (40) 

Self-consistent ( "^-derivable" ) approximations are obtained by selecting a 
class of skeletons in ^[D] and calculating H from Eq. ( p9| ) above, preserving the 
stationarity condition. 

The stationarity of Q[D] has an interesting consequence for the entropy 
S = —d{n/V)/dT. Because of Eq. (|40|), the temperature derivative of the 
spectral density in the dressed propagator cancels out and only the explicit 
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Bose-Einstcin factors need to be differentiated in (|37|), yielding |75|, |7^, [1^, 17 

-Im log D^^{uj, k) 

-Inin(tj,fc)RcD(cj,fc) +5' (41) 



d k dn{u))^ 
(27r)4 dT 
d'^k dn{uj) 
(27r)4 dT 

with 



S' = 



dT 



fAJL^ncnimD^O (42) 



up to terms that are of loop-order 3 or higher. Thus, in contrast to fl, where 
$ contributes already to order in perturbation theory, Eq. (5) with S' = 



is perturbatively correct to order g^. The first two terms in Eq. (41) represent 
essentially the entropy of "independent quasiparticles" , while S' may be viewed 
as the residual interactions among these quasiparticles ||7^ . 

The same simplification holds true in the presence of fermions and, with 
nonzero chemical potential, extends to the fermion density Af — —d{^/V)/diJ, 

. In a self-consistent two-loop approximation one thus has the remarkably 
simple formulae, now for general theories 

S = -ti-J ^^^[lm\ogD-\oj,k)-lmU{oj,k)ReD{uj,k)] 

-2tr J [inilog^-i(^,fc)-ImE(^,fc)Re5(L.,fc)] , (43) 

= ^2trJ-0^^^[lmlogS'\cu,k)^IraJ:{oj,k)ReSico,k)], (44) 

where n(w) = {e^'^-iy-, J{uj) = (e'^t""'^' + 1)"^ and "tr" refers to aU discrete 
labels, including spin, color and fiavor when applicable. 

In gauge theories, the above expressions have to be augmented by Faddeev- 
Popov ghost contributions which enter like bosonic fields but with opposite over- 
all sign, unless a gauge is used where the ghosts do not propagate such as in axial 
gauges. But because ^-derivable approximations do not generally respect gauge 
invariancej^ the self-consistent two-loop approximation will not be gauge-fixing 
independent. It is in fact not even clear that the corresponding gap equations 
( ^9[ ) have solutions at all or that one can renormalize these (nonperturbative) 
equations. 

For this reason, we shall construct approximately self- consistent solutions 
which are gauge invariant and which maintain equivalence with conventional 
perturbation theory up to and including order g'^ , the maximum (perturbative) 
accuracy of the two-loop approximation for For these approximations it will 
be sufficient to keep only the two transverse structure functions of the gluon 
propagator and to neglect ghosts. 

^For this, one would have to treat vertices on an equal footing v^ith self-onorgies, which is 
in principle possible using the formalism developed in Ref. |p5||. 
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For soft momenta, we know that the leading order contribution is given by 
the HTL's, and indeed there is no HTL ghost self-energy. For hard momenta, 
one can identify the contributions to ( ^ ) below order as those linear in the 
self-energies, 

5'-'-^ = So + 2N, f ^ Re Im ^ 



ANNf 



^ J (27r)4 dT ' w2 - fc2 
/ TTT-TT 1 Re S+Im - Re I]_Im [ (45) 



considering now a gauge theory with Ng gluons and Nf fermion flavors. Because 
the imaginary parts of the free propagators restrict their contribution to the 
light-cone, only the light-cone projections of the self-energies enter. At order 
this is exactly given by the HTL results, without having to assume soft uj, k 
[3^ (as in the above example of scalar electrodynamics, see Eqs. (p^), (p^) 

ReU[^\oj^ = fc2) = tltiuj^ = k^) = = ml^, (46) 

2fc Re E^t^^ (w = ±fc) = 2fc S± (w = ±fc) = 2M^ = Af^ , (47) 

and without contributions from the other components of 11^^ and the Faddeev- 
Popov self-energy. 

There is no contribution cx g^ from soft momenta in (^3|) and ( ^ ) so that 
one is left with remarkably simple general formulae for the leading-order in- 
teraction contributions to the thermodynamic potentials expressed through the 
asymptotic thermal masses of the bosonic and fermionic quasiparticles: 

Lb f ) F 

Here the sums run over all the bosonic {B) and fermionic (F) degrees of freedom 
(e.g. 4 for each Dirac fermion), which are allowed to have different asymptotic 
masses and, in the case of fermions, different chemical potentials. 

Turning now to the next order, g^, let us first recapitulate how this ususally 
arises in the thermodynamic potential. Since this is a static quantity, we 
can use the imaginary-time formalism and concentrate on the zero-modes as 



in Sect. 4.2.1 . For uj — the leading contributions to the self-energies are 
fl^(0,A:) = ml, cx {gTf and llt{0,k) = 0. This gives the "plasmon-effect"| 
term 



p(3) ^ -ATgT 



d^k 



D \ '"-D 



In deriving the self-consistent expressions for entropy and density, Eqs. ( p3| ) 
and (Q), we can no longer use this argument to extract the order g^ term, 
for we have first rewritten the thermodynamic potential in terms of real-time 

^Obviously a misnomer. It is caused exclusively by the Debye screening mass, not the 
plasmon mass = plasma frequency ^ mj^. 
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propagators and self-energies and then used stationarity to drop T and ^ deriva- 
tives of the spectral densities hidden in the full propagators. In fact, from the 
second line of (|37|), one can still isolate the zero- mode contribution (^) by 
/ [...] = [.. ]{w = 0), but in ( p3[ ) and we have products of imaginary 

and real parts times statistical distribution functions. 

Indeed, the order 5^ contributions now arise from both soft and hard mo- 
mentum regimes. In S^°^^ such contributions are due to the singular behaviour 
of dn/dT ~ l/w which does not allow us to expand out the self-energy inser- 
tions perturbatively, but on dimensional grounds gives a contribution '-^ li^/^. 
In jS^ard^ ^jjg other hand, where we could expand out the self-energies as 
in Eq. (^5|), contributions arise from NLO order contributions to 11 and E 
themselves. (In the case of the pressure, this did not happen because of the 
stationarity property of the pressure (it).) 

With this insight, we can now formulate an approximately self- consistent 
dressing of the propagators that is in line with the maximum perturbative accu- 
racy of the <I>-derivable two-loop approximation: for soft momenta, we take the 
(gauge-invariant and gauge-independent) HTL expressions for self-energies and 
propagators; for hard momenta, where according to (^) only the light-cone 
limit of the self-energies contribute below order 17^, the correct leading-order 
contribution is still given by the HTL expressions (^), (^7|), but in order to 
include the contributions completely, we also require the NLO corrections to 
the on-light-cone self-energies. The latter can be calculated by standard HTL 
perturbation theory, and the theorems of Ref. [52| ensure their gauge indepen- 
dence. 



5.2.1 HTL approximation 

As a first approximation let us consider one which only uses the HTL expres- 
sions without NLO corrections thereof. We have seen that this gives the correct 
leading-order interaction term c>c g^, some part of the contribution, and in- 
finitely many formally higher-order terms as well, since we are going to use 
( ^3| ) and ( ^^ "non-perturbatively" , i.e. without expanding out in powers of g 
and truncating. We can do so because expressions ( ^3| ) and ( ^ ) are manifestly 
UV finite, for they involve only the derivatives of the statistical distribution 
functions — in (|37|), the T = UV-divergences are contained in the integration 
domain uj —00, where the undifferentiated n's and /'s do not fall off expo- 
nentially. 

Using the HTL expressions in ( p3| ) and considering for simplicity the pure- 
glue case, we obtain two physically distinct contributions. The first corresponds 
to the transverse and longitudinal gluonic quasiparticle poles. 



(50) 



where only the explicit T dependences are to be differentiated, and not those 
implicit in the HTL dispersion laws ijJt(k) and u>£{k). Secondly, there are the 
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Landau-damping contributions which read 

k 



Sk^^ = -iV,|^|d^^{2arg[fc2-c.2 + n,] 



-2Imnt Re [uj^ ~P - fit]-! + arg[P + n^] - ImH^Re [k^ + n^l^^j. (51) 

The usual perturbative g^-contribution (Eq) is contained in the first term of 
Eq. ([50|); all the other terms in Eqs. (|50|) , (pi]) are of order in a small-g 
expansion. 

In the HTL approximation, only the soft plasmon effect ~ g^ contained in 
^soft -g present, which turns out to equal 



.soft _ aP(3) 



(52) 



which is jSs in the case of pure glue. However, this identification requires a 
peculiar sum rule 

A53 ^ 7Vgy'^i|2Imfl*Re(A-i?r^)-Imn,Re(^,-<))| 

= A5f ^ + A5f ^ = , (53) 

which we found to hold numerically by cancellations in more than 8 significant 
digits. Rather unusually (cp. Ref. (57|), this does not hold separately for the 
longitudinal and transverse sector and moreover holds only after carrying out 
both, the frequency and the momentum integrations in (|5^). 

Although the g^ term in iShtl is only a fraction of the full one, it would make 
similar troubles when expanded out perturbatively, throwing away all higher- 
order terms in g: for large enough g, the perturbative approximation would lead 
to an entropy in excess of the Stefan-Boltzmann limit. 5htl, on the other hand, 
is a monotonically decreasing function of jhjj/T ~ g. 

5.2.2 Next-to-leading approximation 

The plasmon term ~ g^ becomes complete only upon inclusion of the next-to- 
leading correction to the asymptotic thermal masses rrioo and Mao- These are 
determined in standard HTL perturbation theory through 

Sml^ik) = ReSnriuj = k) 

f*i /»\ ^^^^ 

= Re 

where thick dashed and wiggly lines with a blob represent HTL propagators for 
longitudinal and transverse polarizations, respectively. Similarly, 

^SMlik) = SJ:+{lu = fc) = Re ( , r'*'^. + , )U=fc . (55) 
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The explicit proof that these contributions indeed restore the correct plasmon 
term is given in Ref. . 

These corrections to the asymptotic thermal masses are, in contrast to the 
latter, nontrivial functions of the momentum, which can be evaluated only nu- 
merically. However, as far as the generation of the plasmon term is concerned, 
these functions contribute in the averaged form 



J dkk n3g(fc)Re 6IIt{(^ = k) 
Jdkkn'j^^ik) 



(cp. Eq. 



and similarly 



^^j2 _ / dkkn'pj^{k)Re2kS^+{Lu = k) 
J dkk (fc) 



(56) 



(57) 



These averaged asymptotic thermal masses turn out to be given by the remark- 



ably simple expressions |17 



^CfTrho, 



(58) 



ZTT zvr 

where C/ = Ng/{2N). Since the integrals in (56) and ( |57| ) are dominated by 
hard momenta, these thermal mass corrections only pertain to hard excitations. 
Indeed, in Sect. |4.3| we have seen that e.g. the plasmon mass at fc = receives 
a different, namely smaller correction, whereas the NLO contribution to the 
Debye mass is even positive and logarithmically enhanced. 

Pending the full evaluation of the NLO corrections to Re (511 and ReJS, it 
has therefore been proposed in Refs. |15, |l6|, |l^ to define a next-to-leading 
approximation through (for gluons) 



Snla = Shtl 



soft 



where 



includes (|58|). To separate soft (fc 



hard 

mo) and hard (fc 



(59) 

2-kT) 



momentum scales, we introduce the intermediate scale A ~ ^/2'KT'mDCA and 
consider a variation of ca = 5 ... 2 as part of our theoretical uncertainty. 

Another crucial issue concerns the definition of the corrected asymptotic 
mass fhao- For the range of coupling constants of interest {g > 1), the correction 
|(5m^| is greater than the LO value m^, leading to tachyonic masses if included 
in a strictly perturbative manner. 

However, this problem is not at all specific to QCD. In the simple g'^^p'^ 
model, one-loop resummed perturbation theory gives 



9'T'{l~-g) 



(60) 



which also turns tachyonic for 5^1. On the other hand, the self-consistent 
solution of the corresponding one-loop gap equation |2^, |30| 



m 



n(m) = l2g^J2 



-1 



K 



(61) 
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Figure 3: Comparison of the lattice data for the entropy of pure-glue SU(3) 
gauge theory of Ref. ||l9| (gray band) with the range of Shtl (solid lines) and 
Snla (dash-dotted lines) for p. = ttT . . . AttT and ca = 1/2 ... 2. 

(properly renormalized) , whose perturbative expansion begins exactly like (|60|), 
is a nionotonic function in g. In fact, it turns out that the first two terms in a 
(m/T)-expansion of this gap equation, 

= g^T^ - -g^Tm, (62) 

TT 

which is perturbatively equivalent to (|60|) , has a solution that is extremely close 
to that of the full gap equation (for MS renormalization scales p, ~ 2ttT) jl^ . 

In QCD, the (non-local) gap equations are way too complicated to be at- 
tacked directly. We instead consider perturbatively equivalent expressions for 
the corrected ifioo which are monotonic functions in g. Besides the solution to 
a quadratic equation analogous to (|6^) we have tried the simplest Pade approx- 
imant = g^T'^ / (1 -I- f 5), which also gives a greatly improved approximation 
to the solution of scalar gap equations. In QCD, our final results do not depend 
too much on whether we use the Pade approximant |l^ or a quadratic gap 
equation Q. 

The main uncertainty rather comes from the choice of the renormalization 
scale which determines the magnitude of the strong coupling constant when 
this is taken as determined by the renormalization group equation (2-loop in 
the following). 

In Fig. 1^, the numerical results for the HTL entropy and the NLA one are 
given as a function of T/T^ with Tc chosen as Tc — 1.14A^jg. The full lines show 
the range of results for Shtl when the renormalization scale /2 is varied from 
ttT to AttT; the dash-dotted lines mark the corresponding results for Snla with 
the additional variation of ca from 1/2 to 2. The dark-gray band are lattice 
data from Ref. Evidently, there is very good agreement for T > 2.5Tc. 
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From the above results for the entropy density, one can recover the ther- 
modynamic pressure by simple integration, P{T) — P{To) = dT'S{T'). The 
integration constant P{To), however, is a strictly nonperturbative input. It can- 
not be fixed by requiring P(T = 0) = 0, as this is in the confinement regime. It 
is also not sufficient to know that limT^oo P — ^'frco by asymptotic freedom. In 
fact, the undetermined integration constant in P{T)/ Pf^cciT) when expressed 
as a function of as{T) corresponds to a term Cexp{— q;7^[4/3|^^ +0(as)]}, 
which vanishes for as — *■ with all derivatives and thus is not fixed by any 
order of perturbation theory. It is, in essence, the nonperturbative bag con- 
stant, which can be added on to standard perturbative results, too. However, 
in P{T) / Pf-ccc{T) this term becomes rapidly unimportant as the temperature 
is increased, as it decays like T"^. Fixing it by P{Tc) — 0, which is a good 
approximation in particular for the pure-glue case because glue balls are rather 
heavy, one finds again good agreement with lattice data for T > 2.5Tc. 

This approach can be generalized ||l6|, [T^ also to nonzero chemical potentials 
Hf, where lattice data are not availableQ Simpler quasiparticle models ^ 
have already been used to extrapolate lattice data to finite chemical potential 
[ |65| . The HTL approach offers a possible refinement, but that has still to be 
worked out. 



6 Conclusion 

Hard thermal loops, the leading-order contributions to self-energies and vertices 
at high temperature and/or density, form a gauge-invariant basis for a system- 
atic perturbative expansion, as long as one does not run into the perturbative 
barrier formed by the inherently nonperturbative sector of self-interacting chro- 
momagnetostatic modes. But in QCD one faces the additional problem that 
corrections to leading-order results are so large for almost all values of the cou- 
pling of interest that they lead to a complete loss of (apparent) convergence. 
However, we have seen that further resummations which keep as much as pos- 
sible of the effects of HTL resummation without expanding in a power series in 
the coupling may lead to results that remain valid down to temperatures a few 
times the deconfinement phase transition temperature. 



^Lattice results exist however for the response of therrnodynarrdc quantities to infinitesimal 
chemical potentials, namely for quark number susceptibilities pil pSl BOl, and the above 
approach has been applied to those, too, by now 
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